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Abstract 

Flavour changing top decays t — ► cZ°, t — > eg and t — > c-) are 
predicted with invisible rates within the standard model and may 
represent a window on new physics. We consider these processes in 
supersymmetric extensions of the standard model and we show that 
observable rates can be obtained only if the SUSY breaking is non 
universal and flavour dependent. 
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1 Introduction 



Flavour changing neutral currents (FCNC) in the Standard Model (SM) are 
absent at tree level and suppressed by the GIM mechanism at one loop. They 
are particularly sensitive to large mass splittings between quarks of different 
generations. 

Within the SM, the absence of tree-level FCNC implies that inside the loop 
only charged currents can mediate the flavour change and therefore the large 
splitting of the third generation in the up sector can only be effective in 
processes with external d-type quarks. 

The smallness of the mixing angles between the first two generations and the 
third one implies that the effects of a heavy top are particularly evident only 
in FCNC processes involving b quarks, such as B — B mixing, radiative B 
decays, etc. 

The extensions of the standard model to multi-Higgs doublet schemes 
allow for the presence of tree level FCNC which can be tested also with 
external up-type quarks and in particular in the FCNC top decays. This 
study was pursued in refs. Jl|, where it was shown that one could achieve 
values of the branching ratios of the order of 10~ 6 , i.e. orders of magnitude 
bigger than the standard model estimates of the order of 10 -11 . 

Tree-level FCNC vertices are also present in supersymmetric extensions 
of the SM. In particular, in the Minimal Supersymmetric Standard Model 
(MSSM), i.e. in the supersymmetric extension of the SM with minimal parti- 
cle content, FCNC vertices involving fermions, sfermions and gauginos arise 
because of the misalignment of fermion and sfermion low-energy mass matri- 
ces. By imposing R parity, these vertices can only contribute to FCNC quark 
decays at one loop level via penguin diagrams in which the intermediate par- 
ticles are gluinosQ The difference with SM W-mediated penguins, and with 
SUSY chargino-mediated penguins, is that gluino-mediated penguins involv- 
ing external up-type quarks are proportional to the mass splitting between 
up-type squarks which closely follows the one of the corresponding non su- 
persymmetric partners. 

The purpose of this paper is to analyse neutral current top decays, 
t — > cZ°, t — > eg and t — > cy, in supersymmetric extensions of the standard 
model. 

We give the estimates for the MSSM with flavour-universal soft SUSY break- 
ing terms at the GUT scale, and also for a more general class of models where 
such a flavour universality does not hold. The universal case was already dis- 
cussed in the literature & Ej| . In this case we have recomputed both charginos 
and gluinos contributions. 

We disagree with the results of refs. || [| mainly in the evaluation of the 
relevant SUSY mixing angles and in a detail of the actual calculation of the 
decay amplitude. 

We find that the SUSY mixing angle between the second and the third gen- 
eration has been over-estimated by at least one order of magnitude. We also 

*We will neglect neutralino contributions, because they are in general suppressed by 
powers of aw/o-s- 
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find a difference in the result for the amplitude which can be traced back to 
the omission in previous papers of the diagrams involving a helicity flip in 
the gluino line, which dominate the branching ratios when the gluino mass 
gets large. 

The signal for the case where the SUSY breaking is universal turns out too 
low to be detectable. 

In the non-universal case, where the soft SUSY breaking is not flavour blind, 
the branching ratio can be as high as 1CT 5 , which is maybe detectable. In 
this case the contributions from intermediate chargino exchange with respect 
to gluino exchange can be generally neglected. In Section 2 we review briefly 
those aspects of the MSSM which are relevant for FCNC. In section 3 we 
present the details of the calculation. The last section is devoted to the 
presentation of the results and to the conclusions. 



2 The MSSM 

The superpotential of the MSSM is given by: 

W = Y£QiH 2 Uj + YiQiH x D) + Y^UE X E] - ^H X E 2 , (1) 

which is the R-parity conserving supersymmetric generalization of the Yukawa 
interactions of quarks and Higgs bosons in the SM. 

The breaking of SUSY is accounted for by introducing at the GUT scale 
Mq the following soft SUSY breaking terms: 



£s(M G ) = gf(m|J^Qi + 4R)3U4(-|J y .4 
+ hi (m|J l{ + t R (mfj I R + m\h\h x + m\h\h 2 

A u 13 Q l Ah 2 + A%Q\d J R h x + A\ 3 L\e R h x + BiMfnha 
+ X - (fhgfCg + m w W T CW + m^B T CB) + h. c. 



+ 



(2) 



where C is the charge conjugation matrix and i,j = 1, 2, 3. 

We first consider a constrained version of the MSSM in which we impose 
universality of the soft breaking terms at the GUT scale: 



(m~ ) = (m~ ) = (to- ) = (to- ) = (m~ ) = m^Sa, 

V Q L Jij V u RJij V dnlij \ L L Jij \ i R ) ij 

^9 = =ffl B = m V2- ( 3 ) 

With these soft breaking terms, the model is defined by six parameters: 

m o, m i/2, A], tan/5, fi and B. (4) 

The request of radiative breaking of the electroweak symmetry reduces the 
parameters of the model to five: 

mo, m i/2, A , tan/5 and sign(/x). (5) 
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The low-energy values of the couplings in eqs. (|l|) and @ are obtained 
by solving the following renormalization group equations 
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(6) 



where t = l/(167r 2 ) log(M,|/Q 2 ) is related to the running scale Q 2 . One has 
also to consider the additional set obtained through the replacement u <-> d. 
In eq. (H) we have neglected the leptonic contribution. We have denoted 
by Mi the gaugino masses and by Cf the quadratic casimir eigenvalues of 
SU(3), SU(2) and U(l) for the representation R, with the further definition 



(7) 



At the electroweak scale, we choose a basis such that the Yukawa cou- 
plings of up-type quarks are flavour-diagonal, 



1 — 1 D> 



KY 



D ■ 



(8) 



where 



rn 



D 



/ \ usin/3 
diag (m u , m c , m t ) = —Y r 



d- 



m 



D 



diag (m d ,m s ,m b ) 



V2 

vcos/3 d 

V2 D ' 



(9) 



v = yv 2 + v\ = 246 GeV and tan/3 = i^/fi- 

In the super-CKM basis, in which the quark-squark-gaugino couplings 
are equal to the quark-quark-gauge boson ones, the squark mass matrices at 
the electroweak scale are given by 



rn- 



2t 2 



m 



m% rrv~ 

d-LL 

2T 2 

rrij m~ 

&LR URR 



(10) 
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where 



M 2 

ml LL = m 2 ~ L + (m u D ) 2 + (3 - 4 sin 2 cos 2/3 

m i M = ~ fiml cot P- V ~^pf AU 

m l RR = m l R + (™>l) 2 + l M2 z sin 2 9 cos 2(3 

2 M 2 

m\ Lh = K ] m 2 ~ L K + (m d D ) - — £ (3 - 2 sin 2 fl) cos 2(3 



n J _ l> cos /3 ,, + . . 



d-LR 



m 



2 2 ' - v5 1 



m 2 - + (m u D Y - -M\ sin 2 6 cos 2/3. (11) 



In the low tan (3 regime, in which the effects of the Yukawa coupling 
of the bottom quark Yj, can be neglected, the effect of Y t in the evolution 
of the soft breaking masses is to induce a mass splitting between squarks; 
however, the mass matrix for up-type squarks remains diagonal ||, On 
the other hand, for large tan (3 Y\, gets large and sizeable off-diagonal terms 
are induced by the evolution in m|. However, it is easy to see from eq. (^) 
that the off-diagonal terms between c and t are proportional to K 2 s, and 
therefore at most of the order of a few percents, to be compared with the 
value of Ti / 6 quoted in refs. 0, |3|]. 

In figure [l] we report for different values of tan (3 the renormalization 
group evolution of the value of the ratio of the off-diagonal Cl — th matrix 
element squared divided by the average diagonal matrix element squared, 
which is an estimate of the mixing angle in the mass eigenstate basis. 

As we shall see in section f|, such a small mixing angle renders the pre- 
diction of the constrained MSSM for rare top decays hard to detect. This 
is essentially because the only FCNC effects in the squark mass matrix are 
induced by the quark mass matrix, and are again protected by a GIM mech- 
anism: they disappear when the down quark masses are degenerate. 

By relaxing the universality constraints one introduces a substantial 
flavour mixing in the theory, resulting in large contributions to FCNC pro- 
cesses [@-[ll|. There are strong constraints on non-universal soft breaking 



terms involving the first two families; however, the off-diagonal squark mass 
terms between c and t are unconstrained by available data on low-energy 
FCNC processes^. This means that we can envisage a situation in which 
there is a large mixing angle between c and t, resulting in large contribu- 
tions to rare top decays, as we shall see in section [|, in which we consider 
generalizations of the MSSM with arbitrary mixing angles in this sector and 
nonuniversal gaugino masses. 



2 In principle, constraints on these mass terms could be obtained from a study of 
chargino contributions to b — > decay in non-universal SUSY J12[ . 
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Figure 1: The renormalization group evolution of the ratio 5 of the off di- 
agonal element cl — t L of the squared mass matrix over the average of the 
diagonal elements of the same matrix in the minimal SUSY model with uni- 
versal coupling as a function ot t — l/(167r 2 ) log(M<|/(5 2 ). The four curves 
refer to different values of tan/5 : 10, 20, 30, 35 (from the upper to the 
lower) . 
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Figure 2: Diagrams for the processes t — > cZ , t — > eg, t — > 07. Dashed lines 
are for squarks. 

3 The calculation 

The calculation was rederived independently. 

The general diagrams contributing to top decays are given in figure ^| for the 
processes t — > cZ°, t — > eg, t — > 07. 

When the flavour changing is a small effect, its contribution can be 



safely estimated in the "mass insertion" approximation [|7J-[TT], where the 
vertices are kept flavour conserving and the flavour off diagonal terms are 
treated perturbatively as mass insertions in the squark propagators. This 
approximation allows to make estimates which can be easily transferred from 
one model to another, but in general cannot be trusted when the perturbation 
is not small. 

This applies in particular to the left-right top squark mixing: its main effect 
is to split considerably the corresponding mass eigenstates and to make the 
transition less affected by the GIM mechanism. 

In general, there are three form factors, current conservation restricts 
them to only two for off shell photons and gluons and one for the on shell case. 
For the on shell Z° case, only two form factors survive. The contribution 
from the third form factor, due to axial current non conservation when quark 
masses are different from zero, vanishes when the Z° is on shell. We have set 
to zero the charm mass. 

With respect to previous calculations, we take into account the effect of a 
large left-right stop mixing which produces a large mass splitting and makes 
the transitions with a gluino helicity flip possible and large when the SUSY 
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breaking scale increases. 

In the appendix we report the full result for the form factors relevant to 
the three different decays for the gluino mediated penguins, which represent 
the dominant contribution when the branching ratios are large. 

4 The results 

In table |1] we report the results for the constrained MSSM for two values of 
tan /3 and one chargino masse near the experimental limit, i.e around 90 GeV. 
Squark and gluino masses are set between 200 and 300 GeV. 
Given the smallness of the mixing angle, the correct result for this case can be 
obtained already in the mass insertion approximation. For moderate values of 
tan P the chargino contribution dominates: the gluino contribution becomes 
comparable only for tan/3 = 35. In any case the effect is not visible in the 
universal SUSY model. 







tan $ = 2 


tan p = 35 






gluino chargino 


gluino chargino 


BR(t - 


* c a) 


10 -i5 10 -n 


1Q -U 1Q -12 


BR(t - 


-> cry) 


10 -17 10 -12 


IO" 13 IO" 13 


BR(t - 




io- 17 io- 12 


IO" 13 IO" 13 



Table 1: The results of the branching ratios for the MSSM with universal 
soft breaking. 



By abandoning the hypothesis of flavour universality of the soft SUSY 
breaking terms, one can explore the result for sizeable values of the mixing 
angles. We have restricted our analysis to the case with a mixing between 
the second and the third generation only. The most general rotation in this 
four dimensional space (scharm and stop, both left and right) is parametrized 
by six angles that we chose as the rotation angles in each of the six possible 
rotation planes. More explicitely, our rotation matrix R is given by: 

R = R 12 R 13 R 1A R 23 R 24 R 3 \ (12) 

where 
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COS 6 12 


sin 9 12 









— sin 9 12 


COS 6*12 
















1 





V 











1 / 



(13) 



and similarly for R , R , R , R and R . We generate uniformly in the 
six dimensional space of rotation angles a mixing matrix in the basis where 
the squark mass matrix is diagonal with definite eigenvalues. Our default 
choice are the values at the lower experimental bounds: 



8 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 

BR x10 



Figure 3: The integral of the events with branching ratio greater than BR 
generated with random rotation matrix for the t — > eg decay with gluino 
penguins for the default choice of squark mass eigenvalues (upper curve) and 
for a choice with a smaller splitting among the eigenstates (lower curve). 



m Sl = 220 GeV, m~ C2 = 260 GeV, 
m~ ci = 90 GeV, m- C2 = 180 GeV, 

where the labels 1 and 2 denote the proper eigenstates of the mass matrix. 

The gluino mass has been set equal to 154 GeV. 

For each choice of the matrix we calculate the branching ratio. 

In figure |3| we plot for the gluon decay with intermediate gluino ex- 
change the percentage of the cases where the resulting branching ratio is 
bigger than a given value. Such a percentage gives an indication about the 
naturalness of the corresponding values of the branching ratio. The maxi- 
mum value of the branching ratio is of the order of 5. 10 -5 and a reasonable 
percentage of the order of 20% is obtained for branching ratios larger than 
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Figure 4: The same as in figure 2 for the t — > cy decay and the default choice 
of the mass eigenvalues. 

3. 10~ 5 , i.e. for a detectable signal. On the same figure is given the case where 
the eigenvalues of stop states have been changed from 90 and 180 to 120 and 
160. The curve becomes steeper because of the more severe cancellations due 
to the GIM mechanism. 

Figures |] and [5] refer to the photon and Z° case, respectively. The 
values are smaller by a factor ranging from 20 to 50 coming basically from 
the different colour charge. 

In figure |6] we report the variation, for a fixed favourable choice of the 
mixing angles, of the branching ratio with the values of the squark or of the 
gluino masses, rescaled by a factor / with respect to their reference values. 
The dependence upon the gluino mass at fixed default values of the squark 
masses appears to be more relevant than the reverse case. 

For each rotation matrix and choice of the mass eigenvalues we can 
reconstruct the original non diagonal mass matrix and in particular we can 
analyse which entries of this matrix should be large to reach the highest 
values of the branching ratio. This can be seen by plotting, for the events 
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Figure 5: The same as in figure 2 for the t — > cZ° decay and the default 
choice of the mass eigenvalues. 



11 



-4 



x 10 
0.3 



t^Cg 



0.25 

BR 

0.2 

0.15 
0.1 



0.05 



o 
o 
o 



8 



10 



12 



14 



Figure 6: The variation of the branching ratio for a particular and favourable 
choice of the mixing matrix with the squark (upper curve) or the gluino (lower 
curve) mass rescaling by a factor /. 
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generated with a random rotation matrix where the branching ratio is higher 
than a given threshold value, the correlation of two off diagonal entries of 
the mass matrix normalized to the average of the diagonal entries. 

In figures ^ we plot the correlation of 5l — t R versus CR — ti for two dif- 
ferent cuts on the branching ratio. Among all possible correlations, only this 
correlation shows a hole in the center (when both values are close to zero) 
which increases with the value of the cut. The condition that at least one of 
the two off diagonal entries c-i — t R or c R — ti, should be large is a necessary- 
one when a large branching ratio is required. This cannot be realized in the 
minimal SUSY model with universal couplings. 

Flavor changing top decays are expected to be visible only in extensions 
of the minimal supersymmetric standard model where the soft breaking is 
flavour dependent and non universal. 

NOTE: While completing this work we have seen a paper by Lopez 
et ai, |hep-ph 9702350| , where the same processes were considered. Their 



conclusions on the yield of FCNC in SUSY models with universal breaking 
are different from ours. 



Appendix 

We report in this appendix the result for the gluino mediated diagrams. 
The effective vertex can be parametrized as follow: 

-i u(p){ P R {F a L q 2 ^ + F b L ^ + G h ia^q v ) + 

P L {F a R q 2 Y + F b R fa» + G R io^q u )} e, u(p + q) (14) 

where P L>R = < - 1 ^ 75 ' > and = e^T a for gluon and e M = e M l for photon or 
Z°. In the photon and the gluon case gauge invariance implies F a = —F b , a 
relation that we have explicitely checked. 
The form factors F a can be written as 

= £ E { R kBRA,t L FUAB) - R\ L , B R A t R F 2 a L (A,B)} 

47r A,B=1 

F R = jSY. { R L B RA,t R FUA B) - R\ R , B R A ,h F 2 a R (A, B)}. 

47r A,B=1 

(15) 

Analogous expressions can be written for F b and G. The matrix R diagonal- 
izes the squared squark mass matrix, and the indices A and B identify the 
eigenvalues and the corresponding eigenvectors: 
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Figure 7: The correlation for the events generated with random rotation 
matrices of the matrix elements of the original off diagonal squared mass 
matrix cl — tn and c R — t L for two different cuts of the BR. 
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In general / can be split in a piece proportional to Cp = N 2 Z X =4/3 and 
one proportional to Cq = N = 3, the latter is present in the gluon case. For 
the Cq terms we report only the expression of the form factor relevant for on 
shell gluon. For the Cp terms we report the result at q 2 7^ which is needed 
for the Z° case. The form factor F b does not contribute to the branching 
ratio for on shell vector bosons. 

The formulas are parametrized in terms of the following couplings of 
the gauge boson V = g, 7, Z° to the left or right up quarks d^X* and of the 
couplings of the gauge boson to the up squarks in the mass eigenstates basis 
a mv {A,B): 



for the gluon d^ 9 = a™ = —g. 



r — ys, 

2 



for the photon dfj 1 = a^f 7 = — g em o . 



and for Z° 



3 



2 cos W 3 
2cos6L 3 



The Z° coupling is different for left and right handed quarks and 
squarks and therefore is not diagonal in the squark mass eigenstates basis, 
a mz {A,B) = [R a mz r] ) BA where 



ciqqz = diag (af Z , a q £ Z , af z , a^ z ) . 



We have used the Feynman rules described in [TJJ]. The results have been 
obtained with algebraic package FORM|l5|. 

For the photon or Z° F*(., .) = F a (., and G..(., .) = (?..(., .)° F 

qqV r / 1 

L(R) °A,B <- F • 

—{B x (m 2 top , m\,m 2 A ) + B (m 2 top , m|, mf )} 



2q 

a qqV (A,B) C F ■ { 
. m top( m2 B ~ m ]) 



(mf op -q 2 ) 2 q 
3 m 2 op (m 2 B - m 2 ) 



{Bo(rn 2 op ,m 2 g,m 2 A ) - B o (0,m^,m%) - (m\ - m 2 B )Co(m 2 A ,m 2 g ,m 2 B )} 



•> 1 — 2 5^-{ B o(m top ,mg,m A ) - B {q ,m B ,m A ) + (m B - m-)C {m A ,m g ,m B )} 

L \ m top 1 I 

1 mjop_ 

l top 



9 7^2 72^2 {^o ( m to P ,m~ g ,m A ) + B x {m top ,m A ,m g )} 

z v'hop 1 n 

1 [m B — Wig) 2 2 2 2 2 2 \ , / 2 2\/~i 1 2 2 2 \l 

ofl 2%i B °( m top' TO 9' to a) - -«o(g , m A , m B ) + (m B - m § )C {m A ,m g , m B )\ 

2 vntop - Q ) 



1 m 2 



Lop 



^{B x (m top ,m A ,mj) - B x (q ,m A ,m B ) - 2(m B - m g )C (m A , m g , m B )} 



2 ( m to P ~ 9 ! 

' 1 rn / 2 2 2 \ 2.^/2 2 2\ 1 i 

,71 2\l B x(Q ,™< A> m B ) - m-C {m A ,m- -,m B ) - -} 
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1 (m B m ) 2 2 2 2 \n 

-» 7 2 2 ^ ^ ^top^g^A) - B (0,m s ,m B ) - (m A - m B )C (m A ,m g ,m B )}} 
z \ m top 1 )Q 

(17) 

a q ff R) S a ,b C f ■ 

- *° p % {B {ml ,m 2 s ,m A ) - B (0, m 2 A , m?)}} 
1 m top TOg _ {B (m 2 op ,m 2 g ,m 2 A ) - B (0,m?,m|) - (m^ - m 2 B )C (m 2 A ,m 2 g ,m 2 B )} 

Z \ m top 1 Jl 

_ ™ to P m 9 {B (m 2 op , m 2 g ,m 2 A ) - B {q 2 , m 2 B ,m A ) + (m\ - m 2 g )C (m 2 A ,ml,m 2 B )} 

\ m top 1 ) 

I m^rn, {Co{m 2 A ^ m%)}} (18) 

z l TO top 1 ) 

It can be verified that the above expression has a smooth limit when q 2 — > 0. 

Gil(AB)^ = G lfl (A,i?) c - = 
a g -gy(A,S) C F • { 

+ m *op7-^ §^{Bo(m 2 top ,ml 7 m 2 A ) - S o (0,m?,m|) - (m^ - m 2 B )C (m 2 A ,m 2 g ,m 2 B )} 

\ m top 1 I 

3 q 2 mtop( m B ~~ TO f ) r D / 2 2 2 \ d/2 2 2 \ , / 2 2\r< I 2 2 2 m 

+ o — 7 — 2 2xT-{ B vy m topi m g > m A) -B {q , m B , m A ) + {m B - m g )C (m A , m- g , m B )\ 

1 \ m top 1 ) 

1 m top q {B x (m 2 op ,m 2 A ,m 2 g ) - B x {q 2 ,m 2 A ,m 2 B ) - 2(m| - m 2 g )C {m 2 A , mj, m%)} 

L \ m top 1 ) 

1 m t0 p {B {m 2 op ,m 2 A ,m 2 g ) + 2B x (m 2 op ,m 2 A ,m 2 g ) - m 2 B C (m 2 A} m 2 g ,m 2 B ) - i}} 

z \ m top 1 I L 

(19) 

G 2L (A,B) C * =G 2R (A,Bf* = 
a qqV (A,B) C F ■ { 

-\ , 2^ 2\ { B Q{rn 2 op , m 2 g , m 2 A ) - B o (0, m 2 g ,m 2 B ) - (m 2 A - m 2 B )C (m 2 A ,m 2 g , m|)} 
z \ m top <7 ) 

2 

- 1 2 m§q 2 \2 { B °( m topi m l' m A) - Bo(q 2 ,m 2 B ,m A ) + (m 2 B - ml)C (m 2 A ,m 2 g ,m 2 B )} 

\ m top 1 ) 

-\ (j m *2^ m A^WB)}} (20) 

z v'hop q ) 

For the gluon G..(., •) = <?..(•> -) Cp + <?..(., .)° G 
Gil(A, B) Cg — Gir(A,B) Cg = 

-g s 5a,b -jC G ■ { 

-m t0 p , B 2 * {B (m 2 top ,m 2 g ,m 2 A ) - B Q (0,m 2 g ,m 2 B ) - (m 2 A - m B )C {m 2 A ,m 2 g ,m 2 B )} 

\ m top) 

+ \^^{Bo{m 2 top ,m 2 gl m 2 A ) + 2B x (m 2 top ,m 2 A ,m 2 g ) - m 2 B C {m 2 A ,m 2 g ,m 2 B ) - ^} 

Z m top Z 

+m top , 2 \2 {Bo(m top ,m s ,m A ) - B (0,m g ,m A )} 

\ m top) 

+ l-^{Bo{m 2 top ,m 2 g7 m 2 A ) + 2B x (m 2 op , m 2 g , m 2 A ) - m 2 g C (m 2 g ,m 2 A ,m 2 g ) - i}} (21) 

Z m top Z 
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G 2L (A,Bf° =G 2 r(A,B) c ° = 

-9s $A,B TpG ■ { 

+ l-^{ B o(mt op ,m 2 g,m 2 A ) -S o (0,m?,m|) - (m 2 A - m 2 B )C (m 2 A ,rn 2 s ,m 2 B )} 

Z m top 

-J^r{ B oK P , m?, m\) - B (0,m 2 s , m 2 A ) + m 2 op C (m 2 s ,m A , m?)}} (22) 

z m top 



Where: 

B (p 2 ,ml,ml) = pole terms + 

— da log (p 2 a 2 + (m 2 — m 2 . — p 2 )a + m 2 .) 
Jo 

B x (p 2 , m^m 2 ,) = pole terms + 

+ / da ct log (p 2 a 2 + (m\ — m\ — p 2 )a + ml) 
Jo 

C {m 1 ,m 2 ,m 3 ) = 



L 



1 da 1 b (m 2 op a 2 + (m 2 - mj - m 2 op )a + m|) 

o (jn 2 p ~ <l 2 ) a + ( m 3 — m i) (q 2 a 2 + (m 2 — m 2 — q 2 )a + m 2 ) 

(23) 

The above integrals contain divergent parts (pole terms) which cancel 
in the finite result because the matrix R is unitary. 

The integrals are performed in the dimensional regularization scheme: in 
particular the reduction of tensor integrals to scalar integrals is performed in 
4 — e dimensions, while the gamma matrices are kept in 4 dimension. The 
final result being finite is indipendent upon the particular regularization that 
we have chosen. 

In terms of the above form factors the expression for the branching 
ratios reads: 

BR = r(tZb + w) with r ^ ~" b + w ^ ~ 1,52 GeV and 



r(t- 


^ c + 7) 


r(t- 


^c + g) 


r(f- 


>c + Z° 



16tt 
1 

16n 
1 m 2 op -m| 

32vr 2ml p 



r 2 R ) |(g2 =0) 



—C F m\ op (G\ + G|) | (g 2 =0) 



{ (Ff + F a R 2 ) \tf=m%) (2< p m| + 2m\ op m\ - Am%) 
+ (FIG L + F^G R ) | (?2=m |) (-12m 3 top m 2 z + 12m top m z ) 
+ (Gl + Gr) \ { <?=m%) (4mt p - 2m 2 top m 2 z - 2m% )} 

(24) 
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